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A STUDY OF THE TRANSIENT BEHAVIOR OF A HIGH SPEED 
TUNNEL-DIODE SWITCHING CIRCUIT USING AN ANALOG COMPUTER 



by 

GEORGE HANNAUER 
Electronic Associates, Inc. 



L ANALYSIS OF SWITCHING CIRCUIT 

The switching circuit is best understood by refer- 
ence to the simplified schematic drawing in Fig- 
ure 1, The voltages V^. and V_ are bias voltages 
provided by an external source, and are assumed 
to be equal in magnitude and opposite in polarity. 
The tunnel diodes TD^ and TD2 are identical. If 
no inputs or loads are connected to the circuit, 
the same current will flow through both diodes. 
If the volt-ampere curve of the diodes always ex- 
hibits a positive resistance characteristic, the 
voltage drops across the diodes will be identical, 
and the voltage at the juncture of the diodes will 
be zero. 
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LOAD 
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Figure 1. Simplified Schematic for Tunnel-Diode 
Switching Circuit 

The characteristic curve of a typical tunnel diode 
is given in Figure 2, Note that the slope is nega- 
tive over the greater portion of the operating 
range (from about 60 to 300 millivolts). The exist- 
ence of this negative resistance region means that 
the equilibrium described in the preceding para- 
graph will be unstable. A very slight disturbance 
will upset the summetry described above, causing 
a very much larger voltage to appear across one 
diode than across the other. This sensitivity to 
small input disturbances enables the circuit to 
function as a high-gain high-speed switch. 



Assume and V_ are initially zero and that V+ 
gradually increases while V_ decreases, such that 
+ V_ = 0. Then the operating point of each 
diode will move from the origin in Figure 2 along 
the characteristic curve toward the peak. The volt- 
age at the junction will remain zero, as before. 
After the peak has been passed, the voltages and 
currents for the two diodes should remain identi- 
cal (by symmetry), and the voltage at the junction 
should remain zero. However, since the diodes 
are now in their negative resistance region, this 
equilibrium will prove extremely unstable. 

Suppose, for instance, that the circuit is unbal- 
anced by a small positive voltage at the input 
terminal. Then tunnel diode #2 will have a slightly 
greater voltage across it than tunnel diode #1, and 
as the bias voltage increases, tunnel diode #2 will 
pass its peak first. Further increases in the volt- 
age across TD2 will result in a decrease in the 
current through it, and hence a decrease in the 
current through TDl as well. Since TDl is still 
in its positive resistance region, a decrease in 
its current will decrease its voltage drop, increas- 
ing the voltage drop across TD2 (since the tunnel 
diodes are acting as a voltage divider, and the 
sum of their voltage drops must equal 2V^). It 
follows that TD2 will move rapidly through its 
negative resistance region, at a speed limited only 
by circuit reactances, and stabilize somewhere 
in the positive resistance region above 300 milli- 
volts, while TDl stabilizes in the positive re- 
sistance region below 60 millivolts. 

Since most of the voltage drop between the two 
bias sources appears across TD2, it follows that 
the voltage at the junction (the output voltage) will 
be positive. A positive input voltage will therefore 
produce a positive output voltage. The output volt- 
age will be of the same order of magnitude as the 
bias voltage V+ , while the input voltage need only 
be large enough to overcome any inherent imbal- 
ance in the diode characteristics. The input volt- 
age can thus be made very small, and the gain of 
the circuit is limited only by the closeness of the 
match between diode characteristics and between 
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Figure 2, Tunnel-Diode Characteristic Curve 



the positive and negative power supplies. If the 
input and output voltages are equal in magnitude, 
as is the case with cascaded logic elements of 
similar design, high gain means that very large 
resistors can be used, resulting in very large 
fan-in and fan-out figures. 

In high-speed logic and arithmetic circuits, the 
voltages and V. will be out-of-phase alternat- 
ing voltages with approximately sinusoidal wave- 
shapes, D-C levels will be superimposed on these 
sinusoids and, in most applications, this bias will 
just equal the amplitude of oscillation assuring 
that the bias voltages do not change sign. The 
alternating bias voltages will then serve as a 
* 'clock," supplying the necessary timing signals 
for high-speed arithmetic and logic. An output 
can be obtained only during that part of the bias 
cycle when the bias voltage is large enough to 
switch one of the diodes into its negative re- 
sistance region. 

The attractiveness of this circuit for digital ap- 
plications lies in the fact that the input need only 



be large enough to unbalance a S5mfimetric circuit. 
If all components were perfectly matched, the 
slightest input disturbance would unbalance the 
circuit in the proper direction, and the theoretical 
gain would be infinite. In practice, the gain will be 
limited by component tolerances. An estimate of 
the gain (or, equivalently, the fan-in/fan-out cap- 
ability) of the circuit can be made from a static 
analysis of the tunnel-diode characteristic and a 
knowledge of circuit tolerances. Such an analysis 
has been made by Chow (4) who predicted large 
gains for closely matched diodes. 

The above analysis is essentially static. For de- 
sign purposes, a djniamic analysis, which takes 
circuit reactances into account, is necessary for 
at least two reasons: 

1, A knowledge of the switching time is neces- 
sary to determine the maximum clock fre- 
quency (i.e., the maximum rate of informa- 
tion transfer) at which the system can be 
operated. 



2, Since the circuit contains two negative-re- 
sistance elements, it may well prove unstable. 
Instead of the switching behavior described 
above, the system may break into oscillations. 

Both the stability and the switching time will depend 
primarily on the inevitable circuit reactances. Fig- 
ure 3 shows the equivalent circuit for the system, 
with circuit reactances and external loading taken 
into account. The principal reactances involved are 
lead inductance and case inductance (which may be 
linnped together) and diode shunt capacitance. Pro- 
vision is also included for coupling the leads to- 
gether with mutual' inductance, which would corre- 
spond to dressing the power supply leads close 
together during construction. A moderate amount 
of coupling proves to have a beneficial effect on 
switching time when the circuit is operated at 
high frequencies. 
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Figure 3. Equivalent Circuit Showing Reactances 

Experimentation with the actual circuit is possi- 
ble, but extremely difficult due to the high fre- 
quencies involved (the range of clock frequencies 
to be investigated is from 100 to 1000 mc.). The 
basic simplicity of the equivalent circuit makes 
it easy to write equations to describe the system 
and thus study the mathematical model (the equa- 
tions) rather than the system itself. 

XL THE MATHEMATICAL MODEL 

The equations describing the system are given be- 
low. These equations are straightforward applica- 
tions of Kirchoff's laws and the known properties 
of resistors, capacitors, and inductors. They may 
be written down immediately upon inspection of 
Figure 3. 

^1 = ^1 ^1-^^ V V'out 



-^2 = h V^D -\ut 



(2) 



(11) 



m 



Equations 1 and 2 are Kirchoff voltage equations 
for the two major loops. Equation 3 is simply the 
statement that the bias or ^*clock" voltages are 
out-of-phase sinusoids with DC levels superim- 
posed. Equations 4 and 5 state the basic volt- 
ampere relation of a capacitor, and 6 and 7 are 
Kirchoff current equations. Equations 8 and 9 state 
the characteristic relation between voltage and 
current in the tunnel diodes (see Figure 2). On the 
computer, these curves will be represented by 
function generators. Equations 10 and 11 are, of 
course, statements of Ohm's Law. 

Equation 11 was omitted for the purposes of this 
simulation, and replaced by the assumption that the 
input current lin was constant. Values from 0.5 
up to 2.0 ma. were tried. The assumption of con- 
stant input current implies a constant-current 
source, which is not very realistic. This simpli- 
fication is not necessary, and equation 11 could be 
programmed directly by adding one additional sum- 
mer and one inverter to the analog circuit. The 
simplification was made in the present study so 
that the results could be compared directly with 
Herzog's digital solution (Ref. 2), since Herzog 
also made this assumption. 

III. THE COMPUTER PROGRAM 

A. Solving for the Derivatives i^^ and 12 

The first step in deriving a computer program is 
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to solve all differential equations for the highest 
derivatives. This has already been done in equa- 
tions 4 and 5, but equations 1 and 2 contain the 
derivatives and and each derivative appears 
in both equations. To solve for and I2, we must 
treat these as a simultaneous system of two equa- 
tions in two unknowns. This system maybe solved 
by determinants, but the resulting computer cir- 
cuit would be unnecessarily complex, A far simpler 
approach is to solve equation 1 for and equation 
2 for giving 

^*2 = r<-v^2V^VV'Vt) 



tional power. This expression, which was obtained 
from a separate digital computer run by curve- 
fitting, can be programmed with standard sub- 
routines. 

On an analog computer, the best approach is 
straightforward use of a standard function genera- 
tor. The diode function generator (DFG) uses biased 
diode networks to approximate a curve by straight- 
line segments. On the TR-48 variable-breakpoint 
DFG, eleven segments are available. The break- 
points (corners) of the curve may be varied so as 
to group many short straight- line segments to- 
gether where the graph curves most sharply and 
use fewer segments where the graph is straighter. 
In this way, an over-all accuracy of about 1% is 
maintained for this particular curve. 



These equations give 1-)^ in terms of ^2 and vice 
versa. Programming them directly leads to an 
algebraic loop (a loop with no integrators), and 
since this is a two-amplifier loop, positive feed- 
back is present. For stability, we must have a 
loop gain of less than one (Ref 5), By inspection 
of the circuit, it is easily seen that the loop gain 
is M^/L^L2. Since this is always less than one 
on physical grounds, the loop is stable and the 
correct values of and I2 may be obtained without 
determinants (Reference 6). 

B. Generating the Characteristic Curve 

The tunnel-diode characteristic in Figure 2 pre- 
sents a problem in function generation no matter 
which simulation method is used. On a digital com- 
puter, the function values can be stored in main 
memory in tabular form, perhaps with an inter- 
polation subroutine to fill in the gaps between tab- 
ulated values. If enough values and/or a sufficiently 
sophisticated interpolation scheme are stored, a 
tremendous demand is made upon the memory. If 
the digital computer is not called upon to solve the 
dynamic equations as well, this is an acceptable 
technique. This was the approach used by Axelrod 
(3) . The simulation was essentially analog, with the 
digital cor ].uter serving only as ^ function gener- 
ator. 



TABLE 1 

SETUP TABLE FOR VARIABLE DFG 



INPUT (VOLTS) 


OUTPUT (VOLTS) 


0.00 


0.00 


0.35 


5.81 


0.85 


9.34 


1.10 


9.93 


1.45 


9.36 


2.85 


3.67 


3. 85 


1.69 


6.20 


0.89 


7.90 


1. 56 


8.55 


3. 20 


9.00 


6.75 


9.50 


12. 00 



Table 1 gives the appropriate breakpoints and the 
corresponding function values. Note that one of the 
DFG's must accept the input -Vqi/SO and produce 
the output +Ij);l/5. Its input-output graph will look 
like Figure 2 reflected around the vertical axis. 
The appropriate table is obtained from Table 1 by 
making all the input voltages negative. The second 
DFG accepts positive inputs and produces negative 
outputs; its graph looks like Figure 2 reflected 
around the horizontal axis. 



Herzog (2> uses a digital computer for the entire 
problem. The characteristic curve is represented 
by an anal5d:ic expression of four terms involving 
a polynomial, two exponential functions and afrac- 



The volt-ampere curve in Figure 2 requires very 
sharp changes of slope. Since these sharp slopes 
exceed the capabilities of the DFG, a potentiometer 
must be placed in the feedback path of the output 
amplifier to increase the gain of the DFG, 
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Figure 4. Analog Computer Circuit Diagram 
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TR-48 POTENTIOMETER ASSIGNMENT SHEET 
P00-P29 



DATE PROBLEM JlJiNNEI^ DlQPE 



POT 


PARAMETER 


SETTING 


STATIC 
CHECK 


SETTING 


NOTES 


POT 


NO 


DESCRIPTION 


STATIC 


OUTPUT 


RUN 


NO 


CHECK 


VOLTAGE 


NUMBER 1 






00 


2Lj5 


.160 


iO.I66 






00 


0 1 




• 160 


-0-I34- 






01 


02 


2M/5 


0.000 


0.00 






02 


03 


aM/5" 


0-000 


0.00 






03 


04 


l/iO/3 


.020 


-0.2O 






04 


05 


l/Ru 


.1X5 


-0.02 






05 


06 


Iin/50 


.020 


+0.2O 






06 


07 


l/lOOCi/9 


.200 


-0.97 






07 


08 


l/lOO C^/9 


.200 


+0.97 






08 


09 


I.C.A07 


.640 


-6.40 




Artificial I.C Test Only 


09 


10 


Ujl/3 


.188 


-1.88 






10 


1 1 




-188 


-1 .20 






1 1 


1 2 




'800 


+0.83 






12 


13 




■800 


-0-67 






13 


1 4 












14 


1 5 


ri/io 


-300 


-1.50 






15 


16 


ri/io 


-300 


+ 1.50 






16 


17 


A/ 500 


•250 


+1.60 


.260 




17 


IB 


Eoc/500 


-260 


-2.60 






18 


1 9 












19 


20 


'/so 1 1 (0) 


.500 


-5.00 






20 


2 i 


'/5o|iiCo)| 


.500 


+5.00 






21 


22 


Vsoo Vo, (o) 


•100 


+1.00 






22 


23 


'/too|Vp2.(o)| 


.100 


-1.00 






23 


24 












24 


25 


Vzo 


.050 


+ 0.49 






25 


26 


Vzo 


.050 


-0.49 






26 


27 












27 


28 












28 


X 


SWITCH # 1 


LEFT 








29 
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Figure 5. Potentiometer Assignment Sheet 



6 



ELECTRONIC ASSOCIATES INC. 
EDUCATION a TRAINING GROUP 
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fc/'/fe3 - ^^23 TliMMEL Diode 



AMP 
NO. 


FB 


OUTPUT 
VARIABLE 


STATIC CHECK 


NOTES 


CALCULATED 


MEASURED 


CHECK PT 


OUTPUT 


CHECK PT 


OUTPUT 


00 




PART OF PFG 










USE INPUT 
RESISTORS V^ITH 04- 


01 


H 


4-Id,/5 












02 




+ 1 


+ U. 


+ 10.00 








03 




10 sin 




/I /A 
-fo.40 








04 


H 


a 

- T , /2OO 




■f 1 -0375 








05 


H 

1 1 






-0-O3IJ 








06 


/ 


-10 COS 


T 1 -ZO 


-lU.UU 








07 


( 
J 


+ 10 sin d/t 


+ 1 .88 










08 




PartofDFG 










US£ INPUT 

pFsisTORS With O? 


09 


w 


-lon /5 




-9.94- 








10 




-I./5 




-5.00 








1 1 




4.T 




+5.00 








1 2 


y 


^ICi / o 




-4.84- 








13 


T 


-Tr. /5 




4-4.84 








14 






+9-68 


-1.00 








1 5 


/ 


+Vd-i /so 


-9.68 


+ LOO 








16 


1 


DO SO 




+ 1.00 








17 


|4 

* 1 


--VoutAo 




-vj. I 0 








18 








-1.00 








19 
















20 


/ 


+ I./5 


-8.30 


4-5.00 








2 1 


/ 


-I2L/5 


+ 6.70 


-5.00 








22 
















23 

















M653 



Figure 6, Amplifier Assignment Sheet 
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TABLE 2 
SUMMARY OF VARIABLES 



Original 
Variable 



1 2t 



Max Value 



500 mv 



50 ma 



50 ma 



2000 ma/ns 



50 ma 



Scale factor 



7^ Volt/mv 
50 

\ Volt/ma 
5 

7 Volt/ma 
5 

Volt/ma/ns 

7 Volt/ma 
5 



Scaled Variable 



[200 J 



/3 = 



5 seconds machine time 
1 nanosecond problem time 



TABLE 3. SUMMARY OF UNITS OF 
MEASUREMENT USED 

Voltage: millivolts 
Current: milliamps 
Resistance: ohms 

Inductance: nanohenries (10""^ henry) 
Capacitance: nanofarads (10~^ farad) 
Time: nanoseconds (10""^ sec) 

Note: The use of the somewhat unusual units for 
inductance and capacitance guarantees that 
the equations I =CV and V = LI will hold in 
this system of units. 



G. Summary of Static Check Conditions 

For static check calculations, convenient values are 
selected for all integrator output signals (in this 
case I]^, l2, Vdi? and VD2) for the driving func- 
tions, V]^, V2, and I^j^, These values may be chosen 
arbitrarily, but they should be small enough to 
avoid amplifier overloads and large enough to pro- 
vide amplifier output signals that can be measured 
with reasonable accuracy. The values chosen were: 

11 = +25 ma. Vj)^ = +50 mv. 

12 = +25 ma. Vd2 = +50 mv. 



R. 



M 



I. 

m 



A 

Edc 
f 



TABLE 4 
SUMMARY OF PARAMETERS 



r = 3 ohms 



8 ohms 



C = 10"^ nanofarads 



L = 0,4 nanohenry 



<VirL^ Try M = 0 for first run. Prob- 
able maximum value around 0. 2 



nanohenry. 



= 1. 0 ma. 



130 millivolts 
130 millivolts 

150 megacycles/second = 0.15 cycles/na- 
nosecond 



w = 27rf = 0.9425 radians/n.s. 
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V]L = +50 mv. V2 = -V-j^ = -50 mv. 



Vl = Lili -MI2 



^in = 



Parameter values are the same as those for the 
first run. The values for all variables may be cal- 



culated from equations 1 
given below. 

=+49.2 ma. 
^Ci = -24.2 ma. 



10. These values are 



ID2 
lC2 



= +49.2 ma. 
= -24.2 ma. 



out 



+ 8 mv. 



The derivatives may be calculated from equations 
1*, 2*, 4, and 5, These are: 

= Vj3^ = -2420 millivolts/nanosecond 

= -207.50 milliamps/nanosecond 

I2 = -167.50 milliamps/nanosecond. 

All variables that have been calculated can now be 
translated into amplifier output voltages. These 
voltages appear on the amplifier assignment sheet. 
They may be checked against values calculated on 
the circuit diagram (to check the programming and 
scaling) and later checked against actual measured 
voltages on the computer (to check the patching 
and the functioning of the components). The check- 
points should also be calculated and measured for 
integrators. (The checkpoint of an integrator is 
minus the sum of its input voltages. It may be read 
out by patching the summing junction of the inte- 
grator temporarily to the summing junction of a 
summing amplifier that is not being used in the 
problem). These calculated values are also tabu- 
lated in the amplifier assignment sheet. 

The initial condition voltages marked **Test" are 
for static test purposes only, and not for the actual 
run. They should be removed prior to the first run. 



-V2 = -^h 

Solving by determinants: 

■'"1 9*2 2 

If we let Vl = 50 mv, V^ = -50 mv, = l.^ =.0.4 
nanohenry and M = 0.2 nanohenry, then = I2 ^ 
+ 250.0 milliamps/nanosecond. 

The output of amplifier 04 should be -11/200 = -1.25 
volts, and the output of amplifier 05 should be +1.25 
volts. 

IV. RESULTS AND CONCLUSIONS 

The switching transients and output wave shapes for 
constant input current are shown in Figures 7 and 
8. Note in particular that at 300 megacycles the 
output voltage goes through an initial oscillation be- 
fore switching. The output is high for only a short 
period of time late in the bias cycle. At 150 mega- 
cycles, the output rises much more sharply and 
gives a good waveshape. 




Figure 7. Output Waveshape of Switching Circuit 
at 150 Megacycle Clock Frequency 



Note that the static test value of the mutual induct- 
ance M is zero. This value was chosen to break 
the algebraic loop in the static test mode, since 
otherwise it would be very hard to trouble shoot. 
However, this static test value does not check the 
algebraic loop itself, and a supplementary test 
should be included with M 7^ 0. For this supple- 
mentary test, we may as well assume that the arti- 
ficial initial conditions on voltages and currents 
are zero, since this part of the circuit has already 
been checked by the main static check. Equations 
1 and 2 then become: 




012345678 
time - nanoseconds 



Figure 8. Output Waveshape at 300 Megacycle 
Clock Frequency 
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Operation at 300 megacycles is possible, but margi- 
nal. Further experimentation with the model indi- 
cates that for reliable operation with sufficiently 
large fan in/fan out capability the circuit should not 
be operated above about 200 or 250 megacycles. 

It is significant that a preliminary pencil-and-paper 
analysis indicated that the system, ought to operate 
at frequencies up to 1000 megacycles. A circuit 
designed to operate at this frequency would not func- 
tion properly and would be extremely hard to 
troubleshoot, due to the circuit loading of the mea- 
suring devices. The computer solution facilitates 
fundamental understanding of the nature of the 
switching process, as well as providing actual so- 
lutions. The ease of parameter changing on the 
analog computer makes it especially easy to in- 
vestigate the trade-off between gain (i.e. fan in/fan 
out capability) and frequency. 

The results agree quite well with Herzog's digital 
simulation, using the Edsac 11 computer at Cam- 
bridge University, and for further analysis of the 
results, the reader is referred to Her zog' s paper 
(2). It is instructive to compare the two simula- 
tions from the point of view of flexibility and con- 
venience. 

a. The digital program requires the equations to 
be v/ritten as a set of first-order equations, solved 
for the highest derivatives, and then converted into 
algebraic equations by the Runge-Kutta or a^ simi- 
lar integration scheme. While standard software 
is available in many installations to simplify this 
task, it is completely eliminated in an analog sim- 
ilation, since the analog solves differential equa- 
tions as differential equations. 

b. The running time for the Edsac 11 was, accord- 
ing to Herzog, about three minutes per solution, as 
opposed to 25 seconds on the TR-48 computer — a 
speed advantage of seven to one. Nothingprecludes 
running even faster than this, since the components 
are operating well within their frequency- response 
limitations. The solution also can be run in repeti- 
tive operation, in which case the solution takes 



only 50 milliseconds ~ effectively zero solution 
time as far as the operator is concerned. 

c. Even more important than running time per se is 
the time required to obtain meaningful and useful 
results. The output of the analog unit is a continuous 
graph, drawn in ink on 11 x 17" paper. . . a form 
of output which 3n engineer can readily observe and 
interpret. By contrast, the digital results were ob- 
tained from an on-line printer, and considerable 
manual effort was required to translate these re- 
sults into graphical form. 

Of course, on-line devices are available for pro- 
ducing graphical readout from a digital computer. 
Some of Herzog's graphical results were obtained 
by oscilloscope photography. Although this method 
eliminates the tedium of point-plotting, it is messy 
and time-consuming, and offers limited resolution. 
Perhaps the best form of readout (certainly the 
most convenient) involves the use of a digital-to- 
analog converter and an analog X- Y plotter. Wliile 
this method is acceptable, it requires expensive 
conversion equipment to translate the data into 
analog voltages for plotting. Such conversion equip- 
ment is unnecessary with the TR-48 computer since 
the signal is an analog voltage in the first place. 

d. Since only about half of the computing capacity 
of the TR-48 computer is used, the simulation can 
easily be expanded to take additional effects into 
account. Herzog's equations assumed a constant 
current input and ignored the fact that the load is 
not purely resistive. These simplifications could 
easily be removed on an analog simulation by 
addition of a few more amplifiers and potentio- 
meters. On a digital computer, any additional 
complexity in the equations would increase the 
running time. 

An attractive alternative, for example, would be to 
simulate an additional tunnel-diode circuit on the 
analog computer and feed the output of the first 
into the second. This would enable the designer to 
determine how flat the output waveshape of the first 
circuit would have to be in order to trigger the 
second successfully. 
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